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ABSTRACT 

Context. Star formation occurs via fragmentation of molecular clouds, which means that the majority of stars born are members of 
binary systems. There is growing evidence that planets might form in circumprimary disks of medium-separation (< 50 AU) binaries. 
The tidal forces caused by the secondary generally act to distort the originally circular circumprimary disk to an eccentric one. Since 
the disk eccentricity might play a major role in planet formation, it is of great importance to understand how it evolves. 
Aims. We investigate disk eccentricity evolution to reveal its dependence on the physical parameters of the binary system and the 
protoplanetary disk. To infer the disk eccentricity from high-resolution near-lR spectroscopy, we calculate the fundamental band 
(4.7 yum) emission lines of the CO molecule emerging from the atmosphere of the eccentric disk. 

Methods. We model circumprimary disk evolution under the gravitational perturbation of the orbiting secondary using a 2D grid- 
based hydrodynamical code, assuming a-type viscosity. The hydrodynamical results are combined with our semianalytical spectral 
code to calculate the CO molecular line profiles. Our thermal disk model is based on the double-layer disk model approximation. We 
assume LTE and canonical dust and gas properties for the circumprimary disk. 

Results. We find that the orbital velocity distribution of the gas parcels differs significantly from the circular Keplerian fashion. The 
line profiles are double-peaked and asymmetric in shape. The magnitude of asymmetry is insensitive to the binary mass ratio, the 
magnitude of viscosity (a), and the disk mass. In contrast, the disk eccentricity, thus the magnitude of the line profile asymmetry, is 
influenced significantly by the binary eccentricity and the disk geometrical thickness. 

Conclusions. We demonstrate that the disk eccentricity profile in the planet-forming region can be determined by fitting the high- 
resolution CO line profile asymmetry using a simple 2D spectral model that accounts for the velocity distortions caused by the disk 
eccentricity. Thus, with our novel approach the disk eccentricity can be inferred from high-resolution near-lR spectroscopy data 
acquired prior to the era of high angular resolution optical (ELT) or radio (ALMA, E-VLA) direct-imaging. By determining the disk 
eccentricity in medium- separation young binaries, we might be able to constrain the planet formation theories. 

Key words. Accretion, accretion disks - Line: profiles - Stars: pre-main sequence - (Stars:) binaries: general - Methods: numerical - 
Techniques: spectroscopic 



1. Introduction 

Star formation occurs via fragmentation of molecular clouds 
causing about 60% of stars to be born as a member of a bi- 
nary system ( Duquennoy & Mayor||199T |. Large initial spe- 
cific angular momentum results in circumbinary disk formation 
around the protobinary, while circumstellar disks (a circumpii- 
mary and a circumsecondary) are formed around the protostars 
for lower initial specific angular momentum (Bate & Bonnell 
|1997[ l. [Bonavita & Desidera| ( |2007[ l demonstrated that the over- 
all frequency of giant planets in binaiies and single stars does 
not statistically differ among planets discovered by radial veloc- 
ity surveys. On the basis of a comprehensive survey for com- 
panions of 454 nearby Sun-like stars, [Raghavan etaL (2010 1 re- 
vealed that both single and multiple stars are equally likely to 
harbor planets. 

Most of the planet-bearing binaries have large separations 
where the planet formation might be unaltered by the com- 
panion's gravitational perturbation. However, stellar multiplic- 
ity might play a major role in planet formation in medium- 
separation (< 50 AU) binary systems. Based on Doppler surveys. 



Eggenberger & Udry (20T0| showed that about 17% of circum- 



stellar exoplanets are associated with binaries. Among this, five 



circumstellar exoplanets are known to date in < 50 AU sepa- 


ration binary system (Queloz et al. 2000 


Hatzes et al. 


2003 


Zucker et al.||2004[ Lagrange et al.||2006 


Chauvin et al. 


2006 


Correia et al.|2008|l. Two exceptional cases are also known, HW 



Virginis ( |Lee etal.| 2009 ) and CM Draconis (Deeg et al. 2008 ), in 
which planets have been detected in circumbinary orbits. Thus, 
planet formation theories, such as core-accretion ( |Bodenheimer| 
'& Pollac k| 1986[ [Pollack et al.|1996) l or gravitational instability 
( Boss 200 l| l, must be able to explain the formation of planets in 
both circumprimary and circumbinary disk environments. 

The circumprimary disk is tidally truncated at 0.35 - 0.5 
times the binary separation, depending on the binary mass ra- 
tio, binary eccentricity, and magnitude of the disk viscosity 
(Artymowicz & Lubow 1994[ l. Owing to the angular momen- 
tum transfer between the disk and the companion (Papaloizou & 
Piingle 1977 ), the disk is truncated, resulting in a greatly reduced 
disk lifetime being available for planets to form in < 50 AU bina- 
ries ( jCieza et al. 2009 ) that is a severe problem for core-accretion 
scenario. In contrast, the formation of gas giant planets by the 
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relatively rapid-action of the gravitational instability might be in- 
duced by the secondary-generated shock waves if the gas cooling 
time is short. Nevertheless, the disk viscosity can heat the disk 
sufficiently to suppress the formation of clumps ( |Nelson|2000| l, 
but with small viscosity, gravitationally unstable clumps can still 
form ( |Bossl2006l l. 

The disk gas feels the companion's periodic perturbation 
leading to strong interaction at the location of the Lindblad res- 
onances. Waves launched at Lindblad resonances carry energy 
and angular momentum from the binary. The disk experiences 
changes in its angular momentum where the waves dampen, re- 
sulting in the development of an eccentric disk (Lu bow|1991a[ l. 
Several mechanism that may cause wave damping have been 
proposed, e.g. shocks that could be effective in colder disks, or 
turbulent disk viscosity acting as a dissipation source and radia- 
tive damping. Both SPH ( Artymowicz & Lubow|1994 | and grid- 
based ( jKley et al.|2008| hydrodynamical simulations have con- 
firmed the eccentricity development in binaries assuming spa- 
tially constant viscosity for the gas as a source of the wave 
damping. Since the orbit of bodies (dust particles or pebbles) 
is perturbed not only by the periodic gravitational potential of 
the binary but the gas drag as well, the development of disk ec- 
centricity might influence the core-accretion processes. 

The maximum size of the building blocks of planetesimals is 
affected by the impact velocity of sub-micron-sized grains in the 
dust coagulation process. According to Zsom et al. (2010), the 
aggregate sizes are lower in eccentric disks than in axisymmetric 
disk environments owing to the increase in the relative velocity 
between the dust particles. We note, however, that an investiga- 
tion of the SED slopes of medium-separation T Tauri binaries by 
|Pascucci et al. ( [2008) showed that the extent of dust processing 
in the disk surface layer and the degree of dust settling in binary 
disks do not differ significantly from those in disks around single 
stars. 

The planetesimal accretion phase, leading to between km- 
sized planetesimals and several lOOkm-sized planetary em- 
bryos, should proceed in an environment where the mutual en- 
counter velocities of planetesimals are on the order of planetes- 
imal surface escape velocities. In this environment, the plan- 
etary embryos could grow quickly by runaway growth mode 
( Wetherill & Stewart|1989 1. Since the runaway growth mode is 
sensitive to the encounter velocities, for increased encounter ve- 
locities, e.g., due to the stirring up of the planetesimal swarm, 
runaway growth might be stopped. [Thebault etal. (2006 ) showed 
that the impact velocity of different-sized planetesimals tends to 
increase owing to the interaction between the companion and 
the gaseous friction in binaries with separations of 10 < abin ^ 
50 AU. Paardekooper et al.l ("2008) found that the planetesimal 
encounter velocities with different sizes could be larger by an 
order of magnitude in eccentric disks than in the axisymmetric 
case. Consequently, planetesimal accretion might be inhibited in 
highly eccentric disks. 

In accretion disks, the presence of the double-peaked emis- 
sion lines are the natural consequence of the gas parcels moving 
in Keplerian orbits around the host star. Huang| ( [T972[ ) and |Smak| 
( |198 1^ presented this in connection with the emission lines of 
Be and cataclysmic variable stars. Horne & Marsh ( 1986 ) inves- 
tigated the emerging line profiles in accretion disks. Since the 
Keplerian angular velocity of gas parcels is highly supersonic in 
accretion disks, the Doppler shift of the line emitted by individ- 
ual gas parcels exceeds the local line-profile width. Summing 
up the line profiles emitted by individual rings of gas parcels, 
and taking into account the radial dependence of the line surface 
brightness, the result is the well-known double-peaked broad 



symmetric line shapes ( |Horne & Marsh] 1986| l. Azimuthal asym- 
metries in disk surface brightness (e.g. density perturbations in 
optically thin lines, or supersonic anisotropic turbulence in satu- 
rated lines), will break the line profile symmetry (Horne 1995jl. 
Gas parcels orbiting non-circularly (i.e., in elliptic orbit) might 
also produce the asymmetric line profiles presented by Foulkes] 
et al.| ( 2004) l for cataclysmic variables, and Regaly et al., (2010) 
for protoplanetary disks. 

In this paper, we investigate the eccentricity evolution of a 
circumprimary disk in a young binary system. We perform an 
extensive parameter study to reveal the dependence of the disk 
eccentricity on several parameters, such as binary and disk ge- 
ometry, and gas viscosity. The hydrodynamical simulations were 
done in 2D by a grid-based parallel hydrodynamic code FARGO 
( Masset|2000| . We model the circumprimary disk evolution un- 
der the gravitational perturbation of the orbiting secondary as- 
suming ff-type viscosity (Shakura & Sunyaev 1973). We calcu- 
late the fundamental band (4.7 fim) ro-vibrational emission lines 
of the molecule ^^C^^O emerging from the disk atmosphere, 
providing a tool to determine the disk eccentricity from high- 
resolution near-IR spectroscopy by means of line profile distor- 
tions. Our thermal disk model is based on the double-layer disk 
model of Chiang & Goldreich ( 1997[ l. Since the velocity distri- 



bution of the gas parcels show supersonic deviations from the 
circular Keplerian one, owing to the eccentric disk state, asym- 
metric molecular line profiles emerge from the optically thin 
disk atmosphere. 

The paper is structured as follows. In the next section, 
we present our hydrodynamical simulation modeling of the 
evolution in general of disk eccentricity. The calculation of 
fundamental-band CO ro-vibrational emission lines emerging 
from an eccentric circumprimary disk are presented in Sect. 3. 
In Sect. 4, we present an extensive parameter study to investi- 
gate the evolution of the disk eccentricity for a wide range of 
binary and disk parameters. Section 5 deals with the comparison 
of our results to other recent simulations, and the observability 
of eccentric signatures. The paper closes with conclusions. 

2. Hydrodynamic disk model 

In our simulations, we use a locally isothermal version of 
FARGO, a publicly available parallel 2D hydrodynamical code 
(Masset 2000). We apply a-type disk viscosity ( Shakura &| 
Sunyaev 1973), assuming a - 0.02 in our model. We adopt 
dimensionless units for which the unit of length is the orbital 
separation of the binary (flbin), and the unit of mass is the mass 
of central star. The unit of time to is obtained from the orbital 
period of the binary, thus to = l/27r, setting the gravitational 
constant G to unity. 

The computational domain is covered by 270 radial and 500 
azimuthal grid cells, and the origin of the grid is on the pri- 
mary star The radial spacing is logarithmic, while the azimuthal 
spacing is equidistant. The disk's inner and outer boundary is at 
O.Olabin and 1 .5flbin, respectively. The initial orbital separation of 
the binary is abin = 1, thus the computational domain contains 
the orbit and the Roche lobe of the system. Open boundary con- 
ditions are assumed at the inner and outer boundaries, i.e the disk 
material is allowed to flow out from the computational domain, 
but no inflow is allowed. The secondary is allowed to accrete the 
material flowing through its Hill sphere with the same rate as 
described in |Kley| ( |1999 ). The secondary feels the gravitational 
pull of the disk, thus its orbital elements can change. 

We use a rotating frame that co-rotates with the binary. The 
binary orbit is initially circular (ebin = 0). The initial surface 
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density profile follows a power-law distribution I.(R) - Sq/? "'^ . 
We use the canonical constant aspect ratio h - 0.05 disk approx- 
imation. The disk self-gravity is neglected because the Toomre 
parameter QiR) - hMt/nR^'L(R) » 1 for our model ( jTbomre| 
[19641 ). 

To model an existing binary system V807Tau with a sole 
circumprimary disk (Hartigan & Kenyon'2003'), the secondary- 
to-primary mass ratio is set to 0.3. The binary eccentricity is ne- 
glected in this particular calculation. The distance unit is taken to 
be 40 AU accordingly to the measured separation ( |Pascucci et al.| 
2008 I. The initial disk mass is set to 2 x 10"^ Mq by the choice 
of S = 17.2 g/cm^ at 1 AU. To prevent numerical instabilities at 
low density regions, an artificial density floor is applied: when- 
ever the disk surface density is below a certain limit, it is reset to 
that limit, which is 10"'^ in dimensionless units. 

Initially, we placed the secondary on a circular orbit into the 
unperturbed disk, in which the surface density had been arti- 
ficially damped. With this damping, the secondary orbits in a 
practically gas-free environment. We use a Gaussian damping, 
i.e., the density is damped as 



Initial state- 



Es + (1 - 2;s)exp 



R - Rn 
o-ti- 



(1) 



where Ss is the surface density at the outer edge of the damped 
region, R^ is the radius where the density damping begins, and 
cTti is the radial extent of the damped region. The damping pa- 
rameters are R^ = 0.6, o-fy - 0.2, and Ss - 10"^. 

2.1. Formation of eccentric disk 

We present our results on the disk evolution through a couple 
of hundred orbits of the binary. The evolution of the 2D surface 
density is shown in Fig.[r| and the azimuthally averaged density 
profile in Fig.|2](panel a). 

After the first couple of binary orbits, double spiral waves 
appear In general, the disk remains approximately axisymmet- 
ric until the ~ 150th binary orbit (Fig.[T]i. Comparing the initial 
and the evolved density profiles in the 50th and 150th orbits, the 
profile becomes less steep, while the density at the disk inner 
edge increases as the disk material is piled up by the secondary 
(Fig.|2] panel a). The density profile is steepened during the sub- 
sequent orbits because of the angular momentum removal by the 
secondary. 

The disk eccentricity is calculated at each individual grid 
cells using the radial and azimuthal velocity components apply- 
ing the equations presented in Regal y et al.| (f2010). To obtain 
the radial disk-eccentricity profiles, the eccentricity distribution 
is azimuthally averaged within the disk radius, which is defined 
as the radius containing 90% of the mass. From the 150th orbit, 
the disk eccentricity rapidly increases (see the evolution of disk 
eccentricity in Fig.[2]panel b). It is expedient to radially average 
the disk-eccentricity profiles at each secondary orbit, in which 
case we can assign a single value to each eccentric state of the 
disk. Since the disk radius and eccentricity profile (Fig.[2] panel 
d) slightly change as the disk and the secondary mutual positions 
vary, their values as a function of time are "noisy". To remove 
this "noise", the radius and eccentricity evolution curves calcu- 
lated at each secondary orbit are smoothed (Fig.|2] panel c). 

At the beginning of the simulation, the disk is truncated by 
the secondary to ~ 15 AU. From the 200th binary orbit, the disk 
begins to expand slightly, and the disk eccentricity starts to grow. 
The disk eccentricity reaches its overall maximum during the 
~ 300th binary orbit, and begins to decline during subsequent 
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Fig. 1. Evolution of the disk surface density shown in six snap- 
shots taken at 0, 50, 150, 300, 450, and 600th orbits of the binary. 
The secondary's counter-clockwise orbit is shown with a white 
dashed circle. 



orbits. The disk reaches a quasi-steady state by the ~ 450th 
binary orbit, in which neither the eccentricity nor the disk ra- 
dius evolves further. By this time, the disk radius is stabilized at 
~ 15 AU. Noteworthy is the correlation of the disk eccentricity 
and radius evolution. 

Regarding the orbital parameters of the secondary, we found 
that they are not changed significantly. Until the 600th orbit, the 
secondary migrates slightly outward to 0.3% of its original dis- 
tance because it gains angular momentum from the disk, while 
the binary eccentricity remains close to zero. 

When the secondary passes close to the outer edge of the 
elliptic disk, a tidal tail develops through which the disk material 
increases the mass of the secondary, but by a negligible amount. 
The evolution of a tidal tail during the 300th orbit is shown in 
Fig.|3] The tidal tail persists for only a quarter of an orbit, and 
reappears in each binary revolution. 

When the disk becomes elliptic, the velocity field departs 
significantly from the circular Keplerian one. The disk is non- 
axysimmetric, not only in density, but in velocity field too. The 
line-of-sight velocity component of the gas parcels calculated 
during the 600th binary orbit is shown in Fig.|4] assuming that 
the disk is seen to be inclined by 60°. Gas parcels that are at 
the same distance to the primary star (shown with white dotted 
circles), hence at the same temperature, have different azimuthal 
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Fig. 2. Panel a: evolution of the azimuthally averaged density profile shown in six epochs calculated at 0, 50, 150, 300, 450, and 
600th binary orbits. The initial density profile is shown with black dotted curve. Panel b: evolution of the azimuthally averaged 
eccentricity shown in five curves calculated at 50, 150, 300, 450, and 600th binary orbits. Only the 0-20 AU region is shown in 
the plot, as the disk is significantly truncated beyond 20 AU. Panel c: evolution of the disk radius (with blue), and the averaged 
disk eccentricity (with black). The curves are smoothed to remove "noise" (see text for explanation). Panel d: variation in the disk 
eccentricity during a full disk precession period. The disk position angle (PA) is measured relative to the position of the orbiting 
secondary. 



velocities. Since the velocity asymmetry at a given radius can be 
larger than the local sound speed (2 - 3kms"') for R < 3 AU, 
the Doppler shift of the lines emitted by a gas parcel exceeds 
the local line width. Consequently, the shape of emission lines 
emerging from the eccentric disk might be distorted ( [Home &] 
Marsh |1986[ l. Therefore, we expect the molecular line profiles 



formed in an eccentric circumprimary disk of a young binary 
system to be asymmetric (see, e.g., |Statler ( 2001| l). 

The elliptic disk precesses retrogradely with respect to the 
secondary's orbit, as if it were a rigid body. The precession pe- 
riod is ~ 6.6 orbital periods of the binary. As a consequence, the 
line profiles are expected to exhibit long-period variations. 

The azimuthally averaged disk eccentricities calculated for 
four different disk position angles (PA) measured with respect 
to the secondary are shown in Fig.|2] panel d. The eccentricity 
of the inner disk (R < 3 AU) shows periodic variation. The disk 
eccentricity reaches its maximum, when the disk is aligned with 
the binary at PA = 180°, i.e., when the disk periastron edge is the 
closest to the secondary. When the disk apastron edge is closest 
to the binary, i.e., for PA = 0°, the disk eccentricity reaches its 
overall minimum. Between these states, when the disk semima- 
jor axis is perpendicular to the binary axis, the inner-disk aver- 



age eccentricity has an intermediate value. As a consequence, 
the line profiles also are expected to display small periodic vari- 
ations on the timescale of the binary orbit. 



3. Fundamental band ro-vibrational CO line profiles 

The spectral lines emerging from an axisymmetric disk with disk 
parcels on circular Keplerian orbits, have symmetric double- 
peaked shapes ( Home & Marsh|1986| ). Analyzing the hydrody- 
namical results, it is reasonable to expect that the optically thin 
atmosphere of an eccentric disk produces non-symmetric molec- 
ular lines as the disk parcels are moving on non-circular orbits. 
We note that several single T Tauri stars display centrally peaked 
symmetric CO line profiles. This might be attributed to an ex- 
tended CO excitation by stellar UV photons (see e.g., [Brittain 
et al. (2007)), in which case the line profiles are also expected to 
be asymmetric. In the following, we do not take into account the 
stellar UV flux, thus our investigation concentrates on systems 
that would produce double-peaked profiles. 

To calculate the strengths of the CO ro-vibrational emis- 
sion lines, we combine the hydrodynamic results with our re- 
cently developed semi-analytical line spectral model (Regaly et 
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Fig. 3. Development of the tidal tail appeared after the secondary 
passes closest to the disk apastron edge during the 300th orbit. 
The secondary's counter-clockwise orbit is shown with a white 
dashed circle. The density enhancement can reach two orders of 
magnitude in the tidal tail. 



al., 2010). In our approach, the thermal disk model is based 
on the double-layer disk approximation of Chiang & Goldreich 
( |1997| l, in which the CO emission lines are formed in the su- 
perheated optically thin disk atmosphere above the cooler op- 
tically thick disk interior The disk inner edge, where the sim- 
ple double-layer assumption cannot be applied, is assumed to be 
a perfectly vertical wall. According to the stellar evolutionary 
tracks ( Siess et al.|2000 l, the luminosity of a 0.3 Mq pre-main- 
sequence star is negligible compared to that of an 1 Mq primary 
(LQ.3Ma/L\Ma - 0-2), hcncc we neglect the irradiation from the 
secondary. We therefore assumed that the disk is heated solely 
by the stellar irradiation of the primary. We note that the accre- 
tion rate measured in our hydrodynamical simulations is below 
2.5 X 10"" Mgyr"', thus viscous heating is negligible compared 
to the stellar irradiation. As a consequence, the disk atmosphere 
is hotter than the disk interior, where a superheated atmosphere 
is formed producing emission molecular spectra. We assume 
that the primary is similar to the one in the system V807 Tau. 
On the basis of the published spectral type of V807 Tau, i.e., 
K7 ( [Hartigan et al. [T994) and assuming that the primary is 
~ 2.5 Myr old, its radius and effective surface temperature are 
taken to be Rt - 1.83 7?o, and r» = 4266 K, respectively, us- 
ing a publicly available pre-main-sequence database (|Siess et al.| 
2000| l. It is assumed that the dust consists of pure silicates with 
a 0.1 /im grain size Praine & Lee|1984| i. The mass absorption 
coefficient of the dust is taken to be /Cy"^' = 2320cm^/g at vi- 
sual wavelengths and ^-^"^^ = 200cm^/g at 4.7yL(m. The dust- 



to-gas and CO-to-gas mass ratios are assumed to be constant 
throughout the disk and being assumed to be = 10"^ and 
Xg -4x 10""*, respectively. 

Taking into account our dust opacity assumptions, the disk 
is optically thick within 3 AU where the CO ro-vibrational fun- 
damental band is excited, regardless of the significant deple- 




Fig. 4. Line-of-sight velocity component of the gas parcels in 
the 600th orbit, assuming a / - 60° disk inclination. The unit 
of contour lines is km The reddish and bluish colors repre- 
sent the receding and approaching part of the disk, respectively. 
Gas parcels at the same distance to the primary (white dashed 
lines) on the approaching side have significantly larger apparent 
velocities than those on the receding side. 



tion that occurs by the time the disk eccentricity has developed 
(see Fig.[2j panel a). The continuum optical depth normal to 
the disk midplane at A.l/um is T4 7^n,(7?) - l,(R)XiiK^''jl^^, which 
drops below 1 (for one side of the disk) if the disk surface den- 
sity is Ed ^ 0.5g/cm^. In our simulation, the surface density 
is S > Ig/cm^ within 3AU (see Fig. [2] panel c), therefore the 
disk interior is optically thick at both 4.7 jjm and optical wave- 
lengths. However, the density drops below this critical value be- 
yond 3 AU (see Fig. |2] panel a), thus the disk is no longer opti- 
cally thick in the outer regions^ Nevertheless, beyond 2-3 AU 
the CO ro-vibrational fundamental band is not excited by the 
primary irradiation, and beyond 5 AU the disk interior has no 
contribution to the 4.7/vm continuum. 

3.1. Asymmetric line profiles 

In advance, here we give a short summary of the physical 
background responsible for the shaping of molecular line pro- 
files emerging in dynamically distorted disk atmospheres. The 
local line profile emerging from the disk atmosphere from a 
given point R,<p in cylindrical coordinates depends on the lo- 
cal Doppler shift. The fundamental line center vo shifts because 
of the apparent motion of the gas parcels along the line-of-sight. 
The shift is significant compared to the line width if the appar- 
ent velocity of the gas parcels exceeds the local sound speedj^ 

' Note, however, that the density can be arbitrarily scaled up result- 
ing in an optically thick disk, as the hydrodynamical simulations are 
independent of the assumed disk mass as long as the change in binary 
orbital elements are small. More details are given in the Discussion. 

- We neglect the effect of turbulent line broadening, thus the intrinsic 
line width is determined by the thermal broadening alone. 
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Fig. 5. V = 1 — > PIO fundamental band CO ro-vibrational line profiles for circular Keplerian (dashed lines) and eccentric 
disk (solid lines) in the 600th binary orbit. The disk extensions were 0.4 AU-40 AU and 0.2 AU-20 AU for panel left and right, 
respectively. The line profiles for inclinations / - 60° and / - 20° are shown with blue and red colors, respectively. The magnitude 
of asymmetry is larger for higher inclination angles. Although the line fluxes normalized to the continuum are clearly different, the 
magnitude of asymmetry remains the same for model #1 and #2. 



Table 1. CO line profile asymmetries for model #1 and #2 



Model 


i 








#1 


20° 
60° 


1.23 
1.10 


1.28 
1.12 


19.6% 
18.2% 


#2 


20° 
60° 


1.33 
1.12 


1.41 
1.15 


21.6% 
22.2% 



" Normalized intensity maximum at the red peak 
Normalized intensity maximum at the blue peak 



Peak-to-peak line profile asymmetry 



Table 2. Flux ratios of models #1 and #2 



^ Disk interior to stellar flux ratio 
Disk edge to stellar flux ratio 
Disk interior -I- disk edge to stellar flux ratio 



|et al.|2007] ), the disk material beyond 3 AU does not contribute 
to the CO line emission in our model. However, material beyond 
this radii out to ~ 5 AU must be considered for the line-over- 
continuum calculation, since the dust in the superheated disk at- 
mosphere contributes considerably to the continuum at 4.7yum. 

It can be clearly seen that the line profiles formed in an ec- 
centric disk (solid lines) are asymmetric, assuming that the semi- 
major axis of the elliptic disk is perpendicular to or at least not 
aligned with the line-of-sight. The line-to-continuum is stronger 
in model #2 than in model #1, owing to the disk continuum (disk 
+ disk edge) of model #1 exceeding that of model #2 (see details 
of the calculated fluxes in Table [2| and the additional hot CO 
emission orbiting between 0.2 AU and 0.4 AU in model #2. The 
depressed disk continuum for model #2 may be due to: (1) the 
rim flux being smaller for a rim closer to the star, since the disk 
rim flux is Fn^ ~ T^^^^rI^^^, where the T^dge ~ R'^'^ (see, e.g.. 

Appendix A of Regaly et al. ( 2010 Q, resulting in ~ R^'^; 
(2) the radial extension of the region close to the rim where the 
disk interior is irradiated directly by the rim (see, e.g. Appendix 
B of Regaly et al. (2010)) is larger for model #1, resulting in a 
stronger continuum for model #1 than model #2. 

The profile asymmetry does not differ significantly between 
the two models. Calculating the peak-to-peak line profile asym- 

To calculate flie Doppler shifts for circularly Keple rian disks and metry as App = 14 - /, |/(0.5[/b + /r] - 1), where and 4 are 

the continuum normalized line fluxes at the red and blue peaks, 
respectively, the asymmetry is App ^ 20% for all models (see 
details in Table[T]). After closer inspection of Fig.|5] the line cen- 
ter is clearly seen to shift toward the peak, which is in excess. 
The magnitude of line center shift also depends on the disk in- 
clination angle and the disk extension. Obviously, the larger the 
inclination angle, the larger the line center shift. Moreover, the 
magnitude of the line center shift is larger in model #2 than in 
model #1 because of the additional emission of hot gas parcels 
on closer orbits in model #2. 



Model 


i 








#1 


20° 


5.13 


2.78 


7.92 




60° 


2.88 


2.78 


5.66 


#2 


20° 


4.05 


1.75 


5.8 




60° 


2.45 


1.75 


4.2 



an eccentric disk, we used Eq. (E.8) and Eq. (4) of Regaly et al. 
( |2010| l, respectively. In the latter case the velocity components 
of the gas parcels (u^(R, (p) and u^(R, (p)) were provided by the 
hydrodynamical simulations. 

Figure |5] shows the line profiles emerging from an unper- 
turbed (dashed line) and an eccentric disk (solid line), assuming 
20° and 60° disk inclination angles. We performed the synthetic 
spectral calculations in two disk models with different disk sizes, 
using the same hydrodynamical output of the 600th binary orbit. 
The computational domain of the synthetic spectral calculation 
was 0.4 AU < < 40 AU and 0.2 AU < < 20 AU for model #1 
and model #2, where the binary separation is 40 AU and 20 AU, 
respectively (left and right panels of Fig.|5]l. In this way, the disk 
edge is at 0.4 AU and 0.2 AU for models #1 and #2, respectively. 
As the CO fundamental band is excited out to 2-3 AU ( |Najita| 



3.2. Formation and variability of line profile asymmetry 

The line profile distortion in eccentric disks occurs owing to the 
same phenomena as for the giant planet-bearing disks presented 
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Fig. 6.V- 1 — >0P10 ro- vibrational CO line-profile variations during a full precession period of eccentric disk calculated in model 
#1 after the 600th binary orbit. For comparison, the line profiles emerging in an unperturbed circularly Keplerian disk are shown 
with dashed lines. Small subfigures show the disk's apparent position with respect to the line of sight. The / = 20° and / - 60° 
inclination angles were assumed, shown in blue and red colors, respectively. The orbits of gas parcels are indicated by green dotted 
ellipses in the subfigures. Panels a and d: show that the asymmetric line profiles emerging in the disk when seen normal to its 
semimajor axis. Panel c: shows that the symmetric line profiles emerging in the disk when seen parallel to its semimajor axis. Panel 
b: shows that the line profiles with a smaller amount of asymmetry emerging from a disk with intermediate position angle. 



in Regaly et al. ( 2010| l. But here the whole disk is eccentric, 
which allows us to give an even simpler explanation of the origin 
of asymmetry. This and the consequences of disk precession are 
explained in this section. 

To provide insight into the cause of the line profile variations, 
we calculated the CO line profiles during a full disk-precession 
period in model #1. Part of the results are shown in Fig.|6] The 
line profile asymmetry is most prominent when the disk semima- 
jor axis is seen perpendicular to the line-of-sight (Fig.|6] panel a 
and d). In contrast, the line profile asymmetry completely van- 
ishes when the disk semimajor axis is seen parallel to the line- 
of-sight (Fig.[6j panel c). 

The formation of the asymmetry can be explained by the 
temperature difference of the gas parcels with the same abso- 
lute values of receding and approaching velocities. In an elliptic 
orbit, the Y velocity component at cylindrical coordinates (R, 0) 
can be given by 



Vy(R,<P) = Vk{R) 



cos(0) + e 
^jl + e cos(0) 



(2) 



where e is the orbital eccentricity of the gas parcels. For simplic- 
ity, we assume that the disk eccentricity is constant with radius. 
Gas parcels at the opposite side of the disk (one at apastron /?ap 
and the other at periastron /?per), have the same absolute value of 
the velocity along the line-of-sight (VyC^ap^^) - -VVC^per, 0)) 
if the equality 



Ran — 



1 



-R. 



(3) 



holds, using Eq. (j2jl. We note that gas parcels at R^^ and Rp^.^ 
are not on the same orbit. Here we align the cylindrical coordi- 
nate system such that at periastron cos(0) = 1, and at apastron 
cos(0) = -1. According to Eq. (l3]l, /?ap < ^pei holds because the 
eccentricity is given by < e < 1 for elliptic orbits. As the gas 
parcels at apastron and periastron are approaching and reced- 
ing, they contribute to the blue and red peaks, respectively. The 
disk is seen by the observer as shown in the subfigure of panel 
a of Fig.|6] Setting the radii of the gas parcels that contribute 
to the red and blue peaks to be /?,ed = R-^p and Ruue = ^per. 
respectively, we find that /?red < ^biue- The temperature of the 
gas parcels at apastron exceeds that at periastron owing to the 
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TaimiR) ~ R ' dependence of the atmospheric temperature on 
distance in the double-layer disk model. Thus, the gas parcels 
receding with speed V make a larger contribution to the red peak 
than the gas parcels approaching with speed -V to the blue peak 
(Fig.|6] panel a). 

Having precessed the disk by 45°, the red peak remains in 
excess, since although the expression Eq. ([3]l does not hold, 
^red < ^biue (Fig.|6j panel b). Although the magnitude of peak 
asymmetry does not change significantly, the departure of the 
line profile shape from the Keplerian one is smaller in the line 
wings. 

By the time the disk has precessed by 90°, the distances of 
the gas parcels receding and approaching with the same speed 
are equal, thus R^g^ = ^biue- Obviously, their contributions to the 
red and blue peaks are just the same, resulting in no asymmetry 
in the line profile (Fig.|6] panel c). 

The disk is seen again perpendicular to its semimajor axis 
in a subsequent 90° of disk precession (Fig. [6] panel d). In this 
phase, the gas parcels at periastron and apastron are now re- 
ceding and approaching, respectively. In this case, /?,ed = Rpei 
and 7?biue = ^ap, thus the distance relation of the receding and 
approaching gas parcels is found to be /?biue < ^red using Eq. 
([3]). The level of red-blue peak asymmetry is nearly the same 
for disks seen at antiparallel position angles perpendicular to 
the semimajor axis (Fig.|6] panel a and d). Although the disk 
eccentricity varies during the disk precession (Fig.|2] panel d), 
this can be explained by its amplitude being small in the regions 
R <2-3 AU where the CO is excited. 

To illustrate the line profile variations during 200 binary or- 
bits, we present trailed spectra of CO in Fig.|7] The periodic 
change in the red-blue peak asymmetry caused by the disk pre- 
cession is clearly visible. In addition, smaller variations in the 
line wings with periods ~ 0.7 times the binary orbital period 
are also visible. This can be explained by small-magnitude disk 
eccentricity variations owing to its dependence on the apparent 
position angle of the secondary with respect to the disk semi- 
major axis. These variations are similar to those that occur in 
superhump binaries, where an additional period is present as a 
beat period between the precessing disk and the binary orbit. 
However, in the superhump case, the period is slightly longer 
than the binary period because of the prograde disk precession 
(Goodchild & Ogilvie 2006). In our case, the period is shorter 
than the binary orbital period because the disk precession is ret- 
rograde. 

To summarize, the line profile becomes asymmetric in the 
quasi-steady eccentric disk state. The period of line profile asym- 
metry is equal to that of the disk precession. Two asymmetric 
phases alternate with each other, and there are two symmetric 
phases amongst them. However, the variation in the red-blue 
peak asymmetry is unlikely to be observed within a decade be- 
cause of the long precession period, which is ~ 250 yr for a 
40 AU separation binary. Nevertheless, as the line wings are also 
subject to variations of period ~ 0.7 times that of the binary pe- 
riod, it might be detectable for ~ 10 AU separation binaries with 
~ 2 Mq primaries. 



4. Parameter study 

To investigate under which conditions the disk will develop ec- 
centricity, we performed an extensive hydrodynamical study. 
The following parameters were varied within physically reason- 
able intervals: the binary mass ratio q, the viscosity parameter 
a, the binary eccentricity ebin, the disk aspect ratio h, the flaring 




-40 -20 20 40 

V (km/s) 



Fig. 7. Trailed spectra ofy = 1 — »0P10 ro- vibrational CO line 
profile at AJ/im emerging from the quasi-steady eccentric disk 
state in model #1 . The line profiles were calculated with cadence 
of 1/20 binary orbit, and the trailed spectra covers 10 binary or- 
bits. Epochs of the line profiles shown in Fig.[6]are marked with 
dashed lines. 

index y, and the disk-to-secondary mass ratio ^disk/sec- In addi- 
tion, we investigated the eff'ects of the open outflow and the rigid 
boundary conditions at the inner boundary of the disk. 

To speed up our calculations, we assumed a larger disk inner 
radius (Ri„ - 0.05flbin) in the parameter study than used previ- 
ously (Rin = O.Olflbin) in Sect. 2. We ran several models with 
O.Olflbin ^ ^in ^ 0.05flbin- The results demonstrated that neither 
the disk eccentricity profile nor the evolution of the average disk 
eccentricity depends significantly on the choice of the disk inner 
radius. The hydrodynamical simulations performed in this pa- 
rameter study cannot be used to calculate the strength of the CO 
emission lines, because when assuming that the distance unit is 
either 40 AU or 20 AU, the hot emitting inner part of the disk is 
not involved in the calculation. This does not prevent the evolu- 
tion of the eccentric quasi-steady state of the disk being able to 
be studied and described in detail. 

4.1. Results for the disk eccentricity and radius 

The radial extension of a circumprimary disk is determined by 
the tidal truncation, whose efficiency depends on the binary mass 
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Fig. 8. Evolution of the disk radii (dashed curves) and the average disk eccentricities (soUd curves) for different models in our 
parameter study: Panel a) models with binary mass ratios 0.05 < q < 0.7, viscosity parameter a = 0.02, and aspect ratio h = 0.05; 
Panel b) models with 0.005 < a < 0.4, q = 0.2, and h - 0.05; Panel c) models with binary orbital eccentricity 0.0 < ebin ^ 0.5, 
q = 0.2, a = 0.02, and h = 0.05; Panel d) models with 0.01 <h< 0.1, q = 0.2, and a = 0.02; flaring disk models with y = 2/7 ai-e 
shown with dot-dashed curves; Panel e) models with disk mass 0.001 < ^disk/sec ^ 0.1, q - 0.2, a = 0.02, and h - 0.05; Panel f) 
models with open outflow and rigid boundary conditions assuming q = 0.2, a - 0.02, and h - 0.05. 



ratio and possibly other physical parameters of the disk itself 
such as the magnitude of the disk viscosity. To characterize the 
time evolution of the disk size, we defined the disk radius as the 
distance measured from the primary star containing 90% of the 
disk material. The disk eccentricity is calculated at each individ- 
ual grid cells as mentioned in Sect. 2.1. To characterize the disk 
eccentric state, the eccentricity distribution was azimuthally and 
within the disk radius radially averaged. Since the 3:1 Lindblad 
resonance point - which is responsible for the eccentricity exci- 
tation ( Lubow|1991a) - is close to the outer disk edge, the disk 
eccentricity is strongly connected to the disk radius evolution. 
Therefore, it is crucial to describe the development of the disk 
radius and the average disk eccentricity in parallel. 

The evolution of the average disk eccentricity tends to follow 
the evolution of the disk radius. The disk is tidally truncated after 



a few tens of binary orbits at a radius R ^ 0.35 - 0.54abin for all 
disk models. Later on, the disk increases in size, then shrinks, 
and finally approaches a size of R ^ 0.35 - 0.4flbin^ when the 
quasi-steady eccentric state has been reached. The disk remains 
circular for a couple of hundred orbits, until at some point the av- 
erage disk eccentricity increases abruptly reaching a temporary 
maximum, which is followed by a decrease. After this jump, the 
eccentricity curves retain their constant values, no further jumps 
or sudden changes occurring within the simulation time. 

The evolution of the disk radius and average disk eccentricity 
differs fundamentally from the above-described general cases for 
thin (h = 0.01) and thick {h - 0.1) models, respectively. While 
in a thick model, the disk does not expand temporarily and the 
average disk eccentricity remains low, for thin models the disk 
eccentricity does not develop at all. Moreover, for significant bi- 



10 



Zs. Regaly et al.: Spectra of eccentric circumprimary disks 



nary eccentricity (ebin ^ 0.2), although the disk expands and 
becomes eccentric temporarily, the average disk eccentricity de- 
cays within ~ 100 subsequent binary orbits. 

The final disk radius is significantly affected by the mass ra- 
tio q of the binary system in models with 0.05 < q < 0.7 and 
a - 0.02. We found that the larger the mass ratio, the smaller 
the final disk radius (Fig.|8] panel a, dashed curves). The final 
disk radius is in the interval 0.35 - 0.45flbini which corresponds 
to 14-18 AU assuming Obin = 40 AU for binary separation. On 
the other hand, the magnitude of the final average disk eccen- 
tricity does not depend on q (Fig. [8] panel a, solid curves). We 
also found that the azimuthally averaged eccentricity profile are 
practically the same in the quasi-steady disk state for all mod- 
els. The timescale of the disk eccentricity evolution shows no 
clear dependence on q. The eccentricity growth rate increases 
monotonically with increasing binary mass ratio in the range of 
0.05 < q < 0.15, while for q > 0.2 the growth rate decreases 
with increasing q. 

To investigate the effect of the magnitude of the viscosity 
on the disk eccentricity evolution, we performed simulations 
with 0.005 < a < 0.04 and q = 0.2. We found that with 
increasing viscosity, the disk reaches a quasi-steady eccentric 
state at larger disk radius. The final disk radii are in the range 
of 0.33abin - 0.42abin, corresponding to 13.5 AU-17 AU (Fig.[8] 
panel b, dashed curves). The magnitude of the viscosity, how- 
ever, has no strong effect on the final value of the average disk 
eccentricity (Fig. [8 panel b, solid curves). Although the final az- 
imuthally averaged disk eccentricity profiles are also very simi- 
lar to each other, a slightly lower value of the final average disk 
eccentricity was found for the a = 0.005 low viscosity model. 
On the other hand, the time required to reach the quasi-steady 
eccentric disk state decreases with increasing values of a. 

To investigate the effect of the initial orbital eccentricity of 
the binary on the disk eccentricity evolution, we performed sim- 
ulations in which the secondary was placed on an initially ec- 
centric orbit with the same apastron distance as used for circular 
models assuming q = 0.2, a = 0.02. The initial eccentricity ebin 
was in the range of 0.05-0.5. We found that the final disk radius 
is also sensitive to the binary eccentricity. As a general tendency, 
the final disk radius for an eccentric binary falls short of that of 
the circular binary case (Fig.|8] panel c, dashed curves). The final 
disk radius decreases with increasing binary eccentricity above 
ebin - 0.2. The disk becomes permanently eccentric only for 
Cbin < 0.2 models, in which the final value of the average disk 



panel c, solid curves). For 
becomes only temporarily 



eccentricity is nearly the same (Fig. [8 
higher binary eccentricity, the disk ' 
eccentric. There is no permanent eccentric disk state at all in 
ebin ^ 0.2 models. The lifetime of the temporary eccentric disk 
state decreases significantly with increasing binary eccentricity 
above eyn - 0.2. 

We modeled circumstellar disks with various aspect ratios in 
the interval 0.01 < /i < 0.1 with q = 0.2 and a = 0.02. We 
found that the choice of the disk aspect ratio substantially af- 
fects the final disk radius (Fig.|8] panel d, dashed curves), which 
increases with increasing h for models 0.03 < h < 0.05, and de- 
creases for thicker disks (0.075 < h < 0.1). For thick (h = 0.1) 
and thin disks (h = 0.01), the disk radii are stabilized early (no 
temporary disk puff-up being present) at ~17 AU and ~ 18 AU, 
respectively. The final average disk eccentricity displays a sim- 
ilar dependence on the aspect ratio as the disk radius (Fig. [8] 
panel d, solid curves). The final value of average disk eccentric- 
ity is at its largest for the h = 0.05 model, but for either larger 
or smaller aspect ratios the disk eccentricity profile is somewhat 
flatter. We note, however, that although the average disk eccen- 



tricity is lower for h > 0.075 than for h - 0.05, the azimuthally 
averaged eccentricity profile of the inner disk (R < 4 AU) is sim- 
ilar. For a very thick disk (h = 0.1), we found that there is no 
temporary high eccentric state, and the disk eccentricity is sig- 
nificantly lower than in the thinner models. We emphasize that 
disk eccentricity does not develop at all in the very thin disks 
{h = 0.01) case. 

We also studied the eccentricity evolution for the flaring disk 
approximation, assuming a y = 2/7 flaring index ( [Chiang &| 
Goldreich|[T997) . The disk aspect ratio (h(R) = h(aun)R'^) was 



h(abm) = 0.03 - 0.1 at the distance of the secondary. In these 
models, the disk radius is stabilized at lower values than in the 
corresponding non-flaring models (Fig.|8] panel d, dot-dashed 
curves). The average disk eccentricity remains below that of 
the appropriate non-flaring ones, except for the thick h = 0.1 
disk case, where the disk eccentricity of the flaring model ex- 
ceeds the appropriate non-flaring one (Fig. [8] panel d, dot-dot- 
dashed curves). Noteworthy is the significantly slower eccentric- 
ity growth rate for flaring models. In general, the azimuthally av- 
eraged disk eccentricity profiles of flared models stay below the 
corresponding non-flaring models, except for the thick {h - 0.1) 
model, where the flared eccentricity profile is above the corre- 
sponding non-flaring one. 

The hydrodynamical equations do not depend on the disk 
mass since So drops out of the hydrodynamical equations ( Kley 
1999). This is true as long as we use the isotherm equation of 
state (p = SCj) and neglect the change in the binary orbital pa- 
rameters. The gravitational potential of the secondary that per- 
turbs the disk, however, is affected by the orbital parameters of 
the secondary, which might be subject to change. The magni- 
tude of the torque exerted by the disk on the secondary indeed 
depends on the disk mass. Thus, the final dynamical state of the 
disk might be influenced by the disk-to-secondary mass ratio 
^disis/sec- To investigate whether the disk mass influences the disk 
eccentricity evolution, we performed several simulations assum- 
ing a 0.0002 - 0.02 disk mass with q = 0.2 and a = 0.02, 
corresponding to 0.001 < g'disk/sec < 0.1. We found that the final 
disk radius does not depend on the initial disk mass (Fig. [8] panel 
e, dashed curves). Similarly, the final value of the average disk 
eccentricity is also insensitive to the disk mass (Fig.|8] panel e, 
solid curves). However, we concluded that the disk eccentricity 
evolution tends to begin at earlier epochs for the more massive 
disks. We note that the azimuthally averaged disk-eccentricity 
profile is practically independent of the disk mass being in the 
above-mentioned range. This can be explained by the changes 
in the orbital parameters of the secondary being negligible, i.e., 
Afl/fl,„ < 0.001 and Aebi„ < 0.001. 

Finally, we investigated whether the choice of the inner 
boundary influences the disk eccentricity evolution. This is an 
important question because a feasible inner boundary condition 
might differ when the disk inner rim is close to or far from the 
primary. There are two extreme boundary conditions regarding 
the outflow mass flux at the inner disk edge: maximum allowed 
and minimum (zero) mass flux occurs for open and rigid-wall 
boundary conditions, respectively. An open boundary condition, 
in which case the mass accretion rate at the inner edge is not 
influenced by the star, seems to be feasible for cases when the 
inner disk is farther away from the star, e.g. when the inner disk 
is cleared by photoevaporation, for example Shu et al. ( 1993[ l 
and |Clarke et al. ( 2001| l. When the disk inner edge is close to the 



stellar surface (i.e. at the magnetospheric radius) the radial ve- 
locity component of the migrating material practically vanishes 
owing to the equilibrium of the stellar magnetospheric pressure 
and the ram pressure of accreting material ( |Koenigl||1991| l. In 
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an extreme situation, the accretion ceases, thus a perfectly rigid 
wall assumption might be valid. Thus, we applied an open out- 
flow (where both the azimuthal and the radial components of the 
velocity are the same in the ghost cells) and rigid inner boundary 
(at which the radial velocity component changes its sign without 
changing its magnitude in the ghost cells) conditions for model 
a = 0.02, q = 0.2. Comparing the results of models with open 
and rigid inner boundary conditions, one can see that both the 
final disk radius and final averaged disk eccentricity profiles are 
independent of the choice of inner boundary condition (Fig. [8] 
panel f, dashed curves). The time needed to reach the quasi- 
steady eccentric disk state is also similar for both boundary con- 
ditions. Noteworthy is the difference in the magnitude of aver- 
aged disk eccentricity value during the temporarily high eccen- 
tricity states. 

5. Discussion 



In his pioneering work, Lubow ( 1991a|b i showed that tidally 
induced eccentricity develops in circumprimary disks of close- 
separation binaries because of the effect of the 3:1 eccentric in- 
ner Lindblad resonance. Kley et al. (2008) numerically investi- 
gated the eccentricity development in circumprimary disks us- 
ing spatially uniform kinematic viscosity. |Lubow| ( |2010[ ) then 
explored the eccentricity growth rate by Unear eccentricity evo- 
lution formulated by Goodchild & Ogilvie ([2006) neglecting 
the effect of viscosity. In our approach, however, the viscosity 
of the disk material is provided by the turbulence within the 



disk (Shakura & Sunyaev 1973 i, where the kinematic viscos 



ity increases with increasing stellar distance. In what follows, 
we compare our results to that of Kley et al. (2008 1 and ,Lubow| 
( |20T0l i. 

5.1. Comparison of simulations 

In line with these studies, our simulations demonstrate that the 
circumprimary disk becomes eccentric with final average ec- 
centricity of e ^ 0.25, assuming a-type viscosity in the range 
0.005 < a < 0.04. In good accordance with the results of Kley 
|et al.| ( p008| l, we found that the disk's final average eccentric- 
ity does not depend on the binary mass ratio. In addition to the 
results of Kley et al. ( 2008[ ), we confirmed these findings for 
mass ratios higher than 0.3. However, the time required to reach 
the quasi-steady eccentric disk state was found to decrease with 
increasing mass ratio only up to ^ - 0.15. Above this mass ra- 
tio, the disk eccentricity evolution decelerated significantly, and 
ever more slowly with increasing q in our simulations. |Kley et al. 
( |200 8 ) found that the increase in the eccentricity growth rate is 
continuous up to <7 = 0.3. This behavior of eccentricity growth 
rate can be explained if tidal forces lower the disk density near 
the resonance, where the eccentricity is excited, and the eccen- 
tricity growth rate then is reduced ( |Lubow|2010[ ). 

Lubow ( |2010) l qualitatively expected that the eccentricity 



growth would accelerate as the viscosity increased because the 
disk is then more likely to spread outwards. Along with the re- 
sults of IKley et al. (2008 1, we confirmed this prediction as the 
disk eccentricity growth rate is found to monotonically increase 
with a. 

Kley et al. ( 2008[ l found that the disk final eccentricity de- 



creases with decreasing values of the dimensionless kinematic 
viscosity up to < 10"^, while below this value the disk eccen- 



viscosity is 2.5 x 10"'' < Vk < 6.5 x lO"*" within the disk radius 
assuming that 0.005 < a < 0.04, and ah- 0.05 non-flaring 
disk. Thus, the values of the kinematic viscosity modeled in our 
simulations are close to the range where the disk eccentricity is 
found to be independent of viscosity by [Kley et al. ( |2008| l. 

The evolution of the disk eccentricity departs significantly 
from the general trends for eccentric binaries with ebin ^ 0.2. 
The 3: 1 Lindblad resonance responsible for the disk eccentricity 
( Lubow|fl991a i Ues farther toward the primary, and the disk is 
truncated at smaller radius in an eccentric binary than in circu- 
lar ones. As a consequence, the excitation of 3:1 Lindblad reso- 
nance might be inefficient, resulting in only a temporary eccen- 
tric state with maximum value of edisk < 0.3 that lasts for only a 
hundred binary orbits. 

We have explored a slightly wider aspect-ratio range than 
Kley et aT] ([2008 ). In their investigation, the disk eccentricity 



tricity is found to be saturated (Fig. 8. of Kley et al. (2008 1). 
In the a-type viscosity approach, the kinematic viscosity is 
Vk = ah^R^^^. Hence, in our models the dimensionless kinematic 



displays a weaker linear increase, the aspect ratio being in the 
range of 0.02-0.06. Lubow (2010 1 concluded that the eccen- 
tricity growth rate increases with increasing sound speed, be- 
ing equivalent to the disk geometrical thickness. We have con- 
firmed this result. However, for very thick disks (h > 0.75) both 
the final disk eccentricity and the growth rate were found to be 
smaller than for the h - 0.05 models, meaning that the aspect 
ratio-eccentricity relation turns above h - 0.75. In addition, the 
disk eccentricity evolution is found to be inhibited for very thin 
disks (h = 0.01). To explain this odd relation between the disk 
geometrical thickness and eccentricity development, further the- 
oretical investigation will be required in which the effect of the 
disk viscosity can also taken into account. 

Noteworthy is the disk eccentricity evolution in flared disk 
models, namely, we have found that a significantly lower aver- 
age disk eccentricity develops in the case of a flaring disk ge- 
ometry than in a flat one. However, there is an exception, i.e., 
the /i(flbin) =0.1 thick flared model, where the resulting disk 
eccentricity is slightly above the corresponding non-flaring one. 
This may occur because a flaring disk is always thinner on aver- 
age than the corresponding non-flaring one in our models. Thus, 
flaring disks are thinner at small R, where they are hence similar 
to the non-flaring models with small h, which also do not have 
high eccentricity. 

5.2. Observability of the line profile asymmetry 

We have shown that the eccentric circumprimary disk of a bi- 
nary, containing considerable amount of gas, produces asymmet- 
ric CO ro-vibrational line profiles because of the eccentric orbits 
{^disk - 0.2) of the CO emitting gas parcels. The largest line- 
profile asymmetry measured between the red and blue peaks of 
the V - 1 — > OP(IO) transition is ~ 20 - 25%, depending on the 
disk inclination angle. However, we keep in mind that the magni- 
tude of the line profile asymmetry depends on the position angle 
of the eccentric disk on the sky with respect to the line-of-sight. 
To observe the largest line profile asymmetry, the disk must be 
aligned perpendicular with its semimajor axis to the line-of-sight 
(Fig.|6j panel a or d). Taking this consideration into account, we 
now present a method to determine the disk eccentricity profile 
by fitting the CO line profile asymmetry. 

Assuming that the disk plane is aligned with the binary or- 
bital plane ( ,Monin et al. 2006| ), which is measurable by the well- 
known astrometric technique ( Atkinson|1966[ l, the disk inclina- 
tion angle can be determined. The disk position angle on the 
sky can also be determined using the spectro-astrometric tech- 
nique (Porter et al. 2005 1, which was applied successfully by 
|Pontoppidan et al.| ( |2008| l. Knowing the disk inclination and the 
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position angle of the eccentric disk on the sky, the model de- 
generacy can be resolved. Taking these considerations together, 
one can find the best-fit disk model in terms of the disk eccentric- 
ity, using our semi-analytical spectral model presented in RegSy] 
et al. ( 2010 1. The Doppler shift of the gas parcels at a given R, ^ 



(model #1) 



point in the eccentric disk seen perpendicular to its semimajor 
axis and inclined by / can then be given by 



Av(R,(f>,i) = 



Vk(R) 



c yi -He(7;)cos(0) 
{sin(0) + (cos(0) + e{R))] sin(0, 



(4) 



where the disk eccentricity profile, e{R), can be approximated by 
a quadratic polynomial function 



ediskW - ao + aiR + 02/?"". 



(5) 



Searching for a best-fit eccentricity profile for model #1, we 
found that oo = -0.0385301, oi = 0.136924, and 02 = 
-0.0092516. Figure|9]shows the line profiles calculated in model 
#1 using a circular Keplerian Doppler shift (solid curves) ap- 
proximating the disk eccentricity using Eq. Q (square symbols) 
befor calculating the Doppler shift. The line profiles emerging 
from the circular Keplerian disk are also shown (dashed curves). 
The line profiles are also calculated assuming / - 20° (blue 
color) and / = 60° (red color). One can see that the models using 
an approximated disk eccentricity (Eq. [5]l perfectly fit the line 
profiles calculated by the models using the hydrodynamical ve- 
locity distributions. Consequently, the disk eccentricity profile 
can be estimated by searching for the best-fit model of an asym- 
metric CO ro-vibrational line profile, with knowledge of both the 
disk position angle and the disk inclination angle. 

A fundamental uncertainty about the observability of eccen- 
tric circumprimary disks whether is there enough time to develop 
the disk eccentricity at all, before the disk is depleted in the 
visous timescale? Since the disk eccentricity evolution ceases 
within ~ 0.2 Myr in all of our models (Fig.[8]l, we expect that 
the eccentric state is occurred well within the disk lifetime. We 
search for the best-fit function of the evolution time required to 
reach the maximal eccentric state as a function of the viscosity 
in the form of f^ax - aa'', finding that a - 0.01 1 and b = -0.47. 
Thus, within limitations of our modeling, the minimum viscos- 
ity required to develop an eccentric disk in a 40 AU separation 
binary is a > 6.8x10 ^2.2xl0-^ and 5x10"^ assuming 1 Myr, 
5 Myr, and 10 Myr gas depletion timescale, respectively. 

5.3. Disk thickness for TTauris 

As we have shown, the circumprimary disk in a binary system 
with ebin ^ 0.2 undergoes a significant increase in eccentricity 
for an aspect ratio in the range of 0.03 < h < 0.075, although 
there might be exceptional cases. The final average disk eccen- 
tricity is found to be significantly lower for thick circumprimary 
disks with aspect ratios of h > 0.1 (Fig. [8] panel d). Moreover, 
for thin disks (h < 0.01) the disk does not become eccentric at 
all. In what follows, we predict the thickness of the protoplane- 
tary disk formed surrounding a T Tauri type star 

We assume that the gas in the circumprimary disk is in ver- 
tical hydrodynamical equilibrium. Assuming that the disk is ge- 
ometrically thin (i.e., h(R) - H/R » 1,7? » R^), the disk aspect 
ratio can be given by 



h{R)^ 



kT{R)R 



(6) 




Fig. 9. Line profiles formed in circular Keplerian disks (dashed 
curves), eccentric disks using the velocity distribution given by 
the hydrodynamical simulations (solid curves), and eccentric 
disk using the best-fit eccentricity profile (squares). The line 
profiles are calculated for disk inclinations / = 20° (blue) and 
/ = 60° (red). 



where M» is the stellar mass, T{R) is the temperature profile 
of the disk, k is the Boltzmann constant, jUg ^ 2.3 is the mean 
molecular weight of the disk gas, and is the proton mass (see 
e.g.pullemond & Dominik] ( |2004 )). In the flaring disk model of 
Chiang & Goldreich \ 1997[ ), the disk interior temperature is 



\ 4 / 



(7) 



assuming hydrostatic and radiative equilibrium. Here the accre- 
tion heating is neglected, i.e., only the stellar irradiation is taken 
into account. In Eq. (j7]i, 6{R) is the grazing angle of stellar irradi- 
ation entering the disk atmosphere, which can be approximated 



2R. ^ ItX'' IR,\-^I' 



(8) 



where Tg = GMtiJ.g/kRt. For T Tauri stars, one can write 
(r,/rg)4/^ ^ 0.007. For an average T Tauri star, if 7? > 0.1 AU, 
6{R) can be approximated by the second term of the right-hand 
side of Eq. ([8]), because this term dominates the equation. After 
some algebra, we get 



h(R) ^ 0.02 



Mo) IaU/ 



(9) 



Figure 10 shows the disk aspect ratio using Eq. (j6]l-(j8]l for three 
different masses of T Tauri type stars. (The stellar parameters 
were taken from |Siess et aT] ( |2000| l assuming 2.5 Myr for the 
stellar age.) As one can see, the disk aspect ratio increases with 
decreasing stellar mass, and is always above 0.01. Although the 
disk temperature is lower for lower mass stars, the squeezing 
effect of the stellar gravitational force is also smaller, result- 
ing in a higher disk aspect ratio. Since the disk aspect ratio is 
~ 0.02 - 0.05 on average, the disks around T Tauri type stars are 
thick enough to be eccentric. We note, however, that the SEDs 
of disks, mapping the dust rather than the gas disk around very 
low mass stars, can be fitted assuming significantly lower aspect 
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Fig. 10. Disk aspect-ratio profiles calculated numerically using 
Eqs. (j6]l-([8]) for three different mass T Tauri type stars. The stellar 
parameters assumed in the calculations are indicated. 



ratios ( Szucs et al.||2010] l. Accordingly, we expect the circum- 
primary disks encircling a T Tauri type young primary star in a 
binary system with ebin < 0.1 to produce asymmetric line pro- 
files owing to the development of disk eccentricity. 

5.4. Limits of our model 

In this study we have used a simple 2D disk model, although our 
assumptions are reasonable to investigate the effect of dynam- 
ically perturbed circumprimary disk on CO ro-vibrational line 
profiles. In the hydrodynamical calculations, we have neglected 
the disk self-gravity. Although the Toomre parameter ("Toomrel 
1 1964 !) is 2 » 1 in our models, the disk self-gravity might 
modify the picture. For global disturbances, such as significant 
disk ellipticity, even though the discs are gravitationally stable, 
pressure and self-gravity can be equally important (Papaloizou 
2002| l. In this case, the parameter measuring the importance of 
self-gravity becomes hQ, which is of the order of unity in our 
simulations. |Marzari et al. (2009) found that the disk self-gravity 
appears to be an important factor, because when included the cir- 
cumprimary disk eccentricity evolution is considerably slower 
than in non self-gravitating disk models. We note, however, that 
Marzari et al. (2009 ) follow the disk eccentricity evolution only 



to ~ 0.01 2 Myr, which is too short to develop significant eccen- 



tricity in our models. In addition, Adams et al. ( 1989] ) found that 



eccentric disk instabilities might also be excited by a process 
called as stimulation by the long-range interaction of Newtonian 
gravity for massive disks. Thus, the inclusion of self-gravity may 
also certainly improve our results. 

For simplicity, we have assumed that the dust and gas tem- 
perature are well coupled in the disk atmosphere. However, 
Glassgold efaL] ( [2004) ) and |Kamp & DuUemond] (j2004 ) found 
that the dust and gas coupling ceases above the disk atmo- 
sphere. As a consequence, the gas in the tenuous region above 
the disk atmosphere might be significantly hotter than predicted 
by the double-layer model. Hence, the slowly rotating distant 
disk parcels above the disk atmosphere could produce a substan- 
tial contribution to the low-velocity part of the line profile, re- 
sulting in a centrally peaked profile. Nevertheless, the centrally 
peaked line profiles are expected to be asymmetric because of 
the supersonic Doppler-shift. To study the effect of the dust-gas 
temperature decoupling on the line profile asymmetry, a more 
sophisticated temperature model is needed. 



We have neglected the irradiation flux coming from the sec 
ondary star. The luminosity of an 2.5 Myr old 0.3 Mq mass sec 
ondary is ~ 20% that of the 1 Mq primary (jSiess et al?]|2000 



The final disk radius is < 0.45flbin on average (Fig.|8]), while the 
disk eccentricity can reach e ^ 0.4 at the disk edge (Fig.|2] panel 
d). Thus, the minimal distance between the disk apastron-edge 
and the secondary is > 0.37flbin, taking into account the eccentric 
disk shape. This distance corresponds to 14.8 AU and 7.5 AU for 
model #1 and model #2, respectively. At these distances, the CO 
ro-vibrational fundamental band is not excited thermally by the 
secondary irradiation. Therefore, the effect of the secondary's ir- 
radiation on the CO line profiles should be investigated only for 
close-separation (< 20 AU) binaries. 

A recent survey of Taurus medium-separation young binaries 
revealed that their secondary infrared fluxes are commensurable 
to their primary ones ( Pascucci et al.|200^ ), hence, young bina- 
ries might have also circumsecondary disks. In some cases, the 
infrared flux of a circumsecondary disk is equal to that of the pri- 
mary, and the common spectra of a system may also contain the 
presumably distorted line profiles of the secondary. To investi- 
gate these systems, more elaborate hydrodynamical simulations, 
incorporating the secondary's disk are needed. 

As mentioned in Sect. 2., the secondary star induces a tidal 
wave during close encounters in each orbit. The gas is signifi- 
cantly compressed inside these temporary waves, as the density 
enhancement can reach two orders of magnitudes, according to 
our hydrodynamical simulations. Since the gas might be heated 
to temperatures high enough to excite the CO, the secondary 
mass flow might have a significant impact on the CO line pro- 
files. Thus, to investigate the effect of the secondary mass flow 
on the line profile asymmetry, we need to incorporate the energy 
conservation in the 2D hydrodynamical simulations, which will 
be the subject of an upcoming paper. 



5.5. Outlook 

Spectroscopy of T Tauri stars detects the emission of molecules 
such as H2O, OH, HCN, C2H2, and CO2. Nevertheless, CO is 
more abundant than any of these molecules by a factor of ~ 10. 
Only H2O could reach the abundance of CO predicted by re- 
cent models that calculate the vertical chemical structure of the 



gas in disk atmosphere (e.g., [Glassgold et al. ( 2004 ), Kamp & 
bullemond (2004), and Woitke et al. (2009)). In some cases, 
such as AATau (C arr & Najita 2008), and both AS205A and 
DRTau ( jSalyk et al. 2008),Tie rotational transitions of H2O 
dominate the mid-infrared (10 - 20 //m) spectra, implying that 
H2O is abundant in disk atmospheres. In contrast to theoreti- 
cal predictions, strong water emission could be the consequence 
of turbulent mixing that carries molecules from disk midplane, 
where they are abundant, to the disk atmosphere ( [Carr & Najita 
2008), or the effects of an enhanced mechanical heating of the 
atmosphere ([Glassgold et al.| 2009 |). While the CO and H2O ro- 
vibrational lines provide information about the inner regions of 
disk out to radii 2-3 AU, their rotational lines are excited in 
the radii rang e 10-100 AU ( [Meijerink et al.[|2008| [Beckwifli &[ 
'Sar gent[[l993 ). As the disk eccentricity beyond 3 AU might be 
large, Cdisk ^ 0.4 (see, e.g., Fig.[2[ panel d), the rotational lines 
are also subject to large distortions owing to the disk eccentric- 
ity. Thus, it is worthwhile to search for asymmetric line pro- 
files emerging from young binaries not only in the ro-vibrational 
spectra of CO and H2O, but also in the rotational spectra of CO 
and H2O. However, as the H2O is heated by stellar X-rays and 
sub-thermally populated beyond 0.3 AU, X-ray heating and non- 
UTE level population treatment is needed to calculate water lines 
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(e.g., |Meijerink et al.| ( |2008l l; |Kamp etaL] ( |20T0l i). The calcula- 
tion of asymmetric mid-infrared rotational molecular line pro- 
files emerging from young binaries will be the subject of a future 
study. Another interesting disk diagnostic tool could be the [OI] 
6300 A line. The stellar UV photons incident on the disk atmo- 
sphere are thought to photodissociate OH and H20, producing 
a non-thermal population of excited neutral oxygen that decays 
radiatively. As a result, [Ol] emission lines emerge from the disk 
atmosphere, as reported by |Acke et al.| ( [2005] l for Herbig Ae/Be 
stars. 



6. Summary 

We have examined the effects of a secondary star on a gas-rich 
protoplanetary disk encircling the primary based on the previous 
work of Lubow|( [T99Talbl l and |Kley etaT| ( i200 8 ). While Kley at 
al. used spatially constant viscosity v, we have assumed a a-type 



viscosity ( |Shakura & Sunyaev 1973 1, which has a spatial depen- 
dence given by V ~ iJ''^. According to our results, the circum- 
primary disk eccentricity begins to increase after the first couple 
of hundreds of binary orbits. We have found that the disk ec- 
centricity grows to a certain maximum value, begins to decline, 
and then stabilizes at e^i^k - 0.2 - 0.35. During the eccentricity 
growth, the disk extends and shrinks back to 0.35-0.45 times the 
binary separation. By the time the disk reaches a quasi-steady 
eccentric state, the surface density profile, the disk truncation 
radius, as well as the eccentricity profile no longer evolve. 

We have found that the final average disk eccentricity is in- 
dependent of the binary mass ratio in the range 0.05 < q < 0.7 
(Fig.|8] panel a), of the magnitude of a in the range 0.005 < 
a < 0.02 (Fig. [8] panel b), of the disk aspect ratio for non-flaring 
models with 0.03 < h < 0.075 (Fig.[t] panel d), of the disk- 
to-secondary mass ratio in the range 0.001 < ^disk/sec < 0.1 
(Fig. [8] panel e), and of the choice of the inner boundary con- 
ditions (open outflow or rigid) (Fig.|8] panel f). In contrast, the 
disk eccentricity is lower for thick disks with aspect ratios above 
h = 0.1, while disk eccentricity does not develop at all for thin 
disks with h < 0.01 (Fig.jSj panel d). The disk average eccen- 
tricity is edisk - 0.2 for a flaring geometry with flaring index 
7 - 2/7, while edisk - 0.35 for non-flaring models assuming 
a reasonable /i(flbin) - 0.05 for the disk aspect ratio (Fig. 10 1. 



We emphasize that the disk becomes eccentric only temporar- 
ily for models assumed to have an eccentric binary orbit with 
ebin ^ 0.2. Thus, the binary's orbital eccentricity protects the 
disk against eccentricity formation that might inhibit the planet 
formation. 

In the eccentric disk, the orbit of gas parcels is non-circular, 
rather than elliptic. The average disk eccentricity is edisk - 0.2 
inside 2-3 AU (see, Fig.|2j panel d), where the ro-vibrational 
fundamental band of CO is excited by the primary's irradia- 
tion. Combining our hydrodynamical simulation with our semi- 
analytic double-layer disk model ( Regaly et al.|2010Jl, we have 
calculated the fundamental band V^l^OPlO ro-vibrational 
line profile of '^C"'0 at 4.75 jum emerging from the super- 
heated disk atmosphere. Since the emission spectra of the CO 
is strongly affected by the supersonic Doppler shift, the CO line 
profiles emerging from the eccentric disk are asymmetric de- 
parting from their well-known symmetric double-peaked form 
( Home & Marsh 1986 1. We have found that the maximum peak- 
to-peak asymmetry is App ^ 20% for our models, which is above 
the detection limit of today's high-resolution mid-lR instruments 
such as CRIRES ( jKaeufl et al. 2004). The peak-to-peak asym- 
metry exhibits periodic variations as the disk precesses on the 
timescale of several times the binary period (Fig.|7]i. As the 



disk precession period is several hundred decades in medium- 
separation binaries, the detection of peak-to-peak asymmetry 
variations is improbable within a decade. On the other hand, 
the slight variations seen on the line wings (Fig. [7]) have periods 
shorter than the binary period, which thus might be detected. 



7. Conclusion 

Our study of both the hydrodynamic and eccentricity evolution 
of the circumprimary disk in medium-separation young binary 
systems has revealed the following findings: 

1 . The quasi-steady eccentric disk state always develops in cir- 
cumprimary disks of young medium-separation (20 -40 AU) 
binaries within the average disk lifetime, if the viscosity is 
between widely accepted values (0.01 < a < 0.1). 

2. The CO line profiles are asymmetric (App ^ 20%) as the 
average disk eccentricity is edisk - 0.2 inside 2-3 AU, where 
the CO is excited by the primary's irradiation. 

3. The orbital eccentricity of binary systems ebin > 0.2 or their 
high/low disk geometrical thickness (h < Q.Ql/h > 0.1) 
might inhibit the development of the quasi-static disk eccen- 
tric state. 

4. The inner {R < 2 - 3 AU) disk eccentricity profile can be 
determined by fitting the observed high-resolution near-IR 
CO line profile asymmetry using a simple 2D spectral model. 

Consequently, taking into account that the eccentricity of pro- 
toplanetary disks might strongly influence planet formation, by 
measuring it we might further constrain the planet formation the- 
ories in medium-separation binaries. 
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